A geometric VOF method for interface resolved phase change and conservative thermal energy advection

We present a novel numerical method to solve the incompressible Navier-Stokes equations for two-phase flows with phase change. Separate phases are tracked using geometric Volume-Of-Fluid (VOF) piecewise linear interface construction (PLIC). Thermal energy advection is treated in conservative form and calculation of VOF fluxes at computational cell boundaries used consistently calculate heat capacity. The boundary as sharp (infinitely thin), which leads discontinuity velocity field across presence difficulty this jump accommodated introduction two-step scheme. has been implemented open source code PARIS validated well-known test cases. These include an evaporating circular droplet microgravity (2D), Stefan problem 3D bubble superheated liquid. accuracy shown results were encouraging. 2D showed excellent prediction volume evolution well preservation its shape. A relative error less than 1% was achieved case, water properties atmospheric conditions. Two cases liquid performed, respective Jacob numbers 0.5 2.15. For final radius bubbles both cases, 7% obtained on coarsest grid, finest.